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Recently quantum tomography has been proposed as a fundamental tool for prototyping a few 
qubit quantum device. It allows the complete reconstruction of the state produced from a given 
input into the device. From this reconstructed density matrix, relevant quantum information quan- 
tities such as the degree of entanglement and entropy can be calculated. Generally orthogonal 
measurements have been discussed for this tomographic reconstruction. In this paper, we extend 
the tomographic reconstruction technique to two new regimes. First we show how non-orthogonal 
measurement allow the reconstruction of the state of the system provided the measurements span 
the Hilbert space. We then detail how quantum state tomography can be performed for multi qudits 
with a specific example illustrating how to achieve this in one and two qutrit systems. 

PACS numbers: 03.67.-a, 42.50.-p 
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I. INTRODUCTION 

With increasing interest in quantum computing, cryp- 
tography, and communication, it is of paramount impor- 
tance that there exist means of benchmarking quantum 
information experiments. A singularly useful tool in this 
regard is Quantum State Tomography (QST), which pro- 
vides a means of fully reconstructing the density matrix 
for a state. The procedure relies on the ability to re- 
produce a large number of identical states and perform 
a series of measurements on complimentary aspects of 
the state within an ensemble. The concept is not new, 
with the first such techniques developed by Stokes |^ to 
determine the polarization state of a light beam. Re- 
cently James et al. Q gave an extensive analysis of 
qubit systems specifically focusing on polarization en- 
tangled qubits, building on earlier experimental work 
but more generally for any number of qubits. We also 
refer the reader to Leonhardt's book Q which gives an 
introduction to some of the concepts and experimental 
techniques of tomography relating to continuous variable 
systems in modern quantum optics. 

It is our aim here to expand on the work of James et 
al. in two ways: firstly, to detail how to perform QST on 
systems of n qudits; secondly, to show how to perform 
QST when access to a full range of single qubit rotations 
and hence the state space is restricted. The first point 
is also motivated with respect to fundamental questions 
regarding non-locality in higher dimensions ||^, ^ as well 
as quantum information processing with improved secu- 
rity for Quantum key Distribution ^ and the need 
to characterize these larger quantum states. The second 
point provides a much larger cross-section of the physics 



community with the possibility of performing QST. 



II. 1 QUBIT 

To start with we will first introduce the Pauli oper- 
ators using the group theoretical definition of them as 
generators. This is not crucial, though facilitates the 
procedure of going to higher dimensions with more sub- 
systems without confusing notation changes. Hence, we 
can write a complete Hermitian operator basis for the 
qubit space: 
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corresponding to the 2x2 Identity operator, Aq, and the 
generators of the SU(2) group, Xj, j — 1,2,3. The reason 
for denoting these with Xj will become apparent as we go 
to higher dimensions. For a single qubit we can always 
write the density matrix as 
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As the generators of SU(2) are all traceless operators, the 
normalization of the density matrix p2 requires rg set to 
one, leaving the other parameters rj=i_...^3 constrained 
only by rf+r'^+r'^ < 1. The terms, rj, can be determined 
from the expectation value of the operators such that 
Tj = (Xj) = Tr[/52Aj]. Thus the single qubit density 
matrix has the form. 
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Theoretically, only three measurements are required to 
define the qubit density matrix. The fourth measure- 
ment, Ao, is practically necessary, as it allows renormal- 
ization of the count statistics to compensate for vari- 
ous experimental biases. The experimental data and the 
calculation of the expectation values (Xj) may lead to 
negative eigenvalues for the density matrix even though 
Tr [P2] = 1. This is due to the intrinsic uncertainty in 
experiments, however the mathematical expression (|^) 
allows such nonphysical states (without the constraint 
r1 + r2 -\- ri < 1 ). By using a maximally likelihood 
technique [gj, a physical density matrix can be derived. 

We note that though the SU(2) generators described 
above do not correspond to any physical state, we can al- 
ways write these operators in conjunction with the Iden- 
tity Ao as a linear combination of physical basis state 
density operators. In spin systems this Pauli group pro- 
vides a perfectly reasonable set of observables, however 
in optics this is not the case. In optics a more common 
example could be the polarization basis. 
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H) = |0), \V) = 
i\l)]/V2. The three 



where, in the computational basis 
|i?) = [|0) + |l)]/V2and \L)^m 
orthogonal measurements are \H),\D) and \L) (depicted 
in Fig (I)). 

Regardless of what orthogonal measurements we choose. 




FIG. 1: Schematic of measurements on the characteristic 
sphere (e.g. Poincare or Bloch ) for qubit quantum state to- 
mography In a), an orthogonal set of |0), |-D) and \L) is shown 
while in b) a non-orthogonal set |0), \9+) and |(^+) is shown. 

we can always write Xj = ajkPk for some other set 
of operators pk- State tomography may then be per- 
formed by measuring the expectation values Ujk — (pk) ~ 
Tr[/52Pfc]. 



A. Non orthogonal state tomography 

In the state tomography that has been previously dis- 
cussed we had assumed that we could measure observ- 
ables at orthogonal points on the characteristic sphere. 
(For instance \H) ,\D) ,\L) in Fig (|]).). In many practical 
situations the method of achieving these measurements 



could be a single qubit rotation followed by a measure- 
ment on |0), more explicitly, single qubit rotation would 
be necessary from |0) + |1) and |0)+i|l) to |0). One could 
envisage many practical situations where it is difficult to 
perform these large single qubit rotations to the |0) state. 
Does this mean that state tomography can not be per- 
formed? The answer is no, state tomography can also 
be performed if one has access only to a small solid an- 
gle on the characteristic sphere. For ideal measurements, 
one still needs to make a set of three measurements that 
project onto |0) and 
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where 9, (j) can be small. Thus we only require a small 
perturbation about some accessible point on the char- 
acteristic sphere (see Fig(^). This observation is likely 
to be important in experiments where qubit rotation is 
more demanding than measurement in the logical basis, 
such as flux qubit systems. 

Naturally, as the measurement axes tend further away 
from orthogonal, the uncertainties for a fixed number 
of measurements will grow accordingly, or alternatively, 
achieving a target uncertainty in the state reconstruction 
will require a larger number of measurements. 

Consider arbitrary states, I'ipv), such that a projection 
measurement is represented by Ay = The count 

statistics arise from a series of these measurements. Cor- 
respondingly the average counts from a series of measure- 
ments will be 



(7) 



where A/" is a constant that will be dependent on exper- 
imental factors such as detection efficiencies. The mea- 
sured counts, Uj/, are statistically independent Poissonian 
random variables and hence we assume that they will sat- 
isfy 



(8) 



This now allows us to consider how these statistics will 
vary with respect to the nonorthogonal measurements. 

The difference in count statistics when measuring with 
orthogonal states and when using nonorthogonal states 
will be proportional to the overlap of the two states . 
We now denote the measurement statistics resulting from 
projecting onto one of a set of nonorthogonal states, \'4''v)i 
as n'jj. Hence we find that the counts for nonorthogonal 
measurements are related to the orthogonal in the fol- 
lowing manner. 



with the errors appropriately scaled and given by 
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The counts and the errors all revert to the orthogonal 
case as ^ 1- 



III. GENERALIZATION TO QUDITS 

We introduced the qubit tomography in terms of the 
SU(2) generators. Let us now consider a state with d 
levels. Firstly we prepare the generators for SU((i) sys- 
tems and thereby construct the density matrices for a 
qudit system. For convenience we use the su algebra but 
we will denote the algebra for a d-dimensional system 
as su(d). The generators of SU(d) group may be con- 
veniently constructed by the elementary matrices of d- 
dimension, {e^|fc,j = 1 . . .d}. The elementary matrices 
are given by 



1 < j^, /i < d, 



(11) 



which are matrices with one matrix element equal to 
unity and all others equal to zero. These matrices satisfy 
the commutation relation: 



This pd is a density matrix of dimension d, a qudit, and 
the coefficient tq is one for the normalization. The con- 
dition Tr [pI] < 1 requires Y.f=i^ <d(d-l)/2. 

Now let us extend these results to n-qudits. It was 
shown that for multiple qubits we only had to consider a 
space of operators defined by the tensor product of the 
generators, SU(2) (g) SU(2) (g)...(g) SU(2) where we have 
included Aq (the normalized identity matrix) with the 
normal SU(2) generators^. For two qudits, a density 
matrix p2 d, which has dimension c?^, can be expanded 
similarly. All combinations of the tensor products of the 
A-matrices (complimented with Aq), Xji ^ Xj2, are lin- 
early independent to each other. Hence, the expression 
of the density matrix p2 a may be written in terms of 
A-matrices, 
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Similarly this expression can be generalized to density 
matrices of n qudits, that is 



[e;-,ef] = 5k,j(^i - Sue']. 
There are d{d — 1) traceless matrices, 



-i(e^ - ei), l<k<j<d 



(12) 

(13) 
(14) 



which are the off-diagonal generators of the SU(d) group. 
We add the d — 1 traceless matrices 



Vr = 



r{r + 1) 



(15) 



as the diagonal generators and obtain a total of d^ — 1 
generators. SU(2) generators are, for instance, given as 
{X = el = el + el Y ^ Pl = -i{e\ - e\),Z 



We now define the A-matrices, this is how we labelled 
the Pauli matrices in Eq.dll), 



Aj2_l 



= e 



(16) 
(17) 
(18) 



which, as shown previously, produce the X, Y, Z opera- 
tors of the SU(2) group and so on for higher dimensions. 
In conjunction with a scaled d-dimensional identity op- 
erator these form a complete hermitian operator basis. 

It is then straightforward from Eq.(||) to see that a 
density matrix pd can be a linear combination of the gen- 
erators as 



3=0 



(19) 



A,„. (21) 



jl...jn=0 

The tomography on such a state is only restricted by the 



{A,jl ® A,j2® ... ® Xjn} 



Particle Dimension (d) and Number of Operations (d^i), ' ji' 



(l,2,...,dl) 



FIG. 2: The measurement scaling for tomography on n-qudits 
results from the necessity to measure every basis state on 
every subsystem in every permutation. The measurements 
scale as d^" — 1 where d is the particle dimension, e.g. d = 2 
for a qubit, and n is the number of particles. 

patience of the experimentalist to determine the expec- 
tation values for the systems observables. 



rji...jn — (A 



Xjn) 



(22) 



There will we d^" — 1 measurements required if we assume 
perfect detection. Fig.(||) illustrates the scaling catastro- 
phe that occurs for multiple parties of higher-dimensional 
states. The key concept in both the extension to higher 
dimensional states and to more subsystems is that for 
each subsystem we need to measure every basis state on 
every subsystem in every permutation. 

However if some structure is known about the state, 
then the number of measurements can be reduced. For 
example if we are confident that we are only ever deal- 
ing with a pure state then the number of measurements 
is significantly reduced and the scaling of measurements 
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more so. QST for two qubits normally requires 15 mea- 
surements. If we know this state is pure this is reduced 
to 6: 3 on the diagonal; and 3 on the anti-diagonal. (In 
the case where we know the state to be, say, one of the 
Bell states, then this is reduced further to just 2). So 
in general for pure states we only require 2(d" — 1) mea- 
surements to reconstruct the density matrix. 

The principle of nonorthogonal state tomography car- 
ries through to the higher dimensional cases in exactly 
the same way that it does for normal tomography using 
orthogonal states as do the considerations with respect to 
errors. Also a detailed discussion regarding the sources 
of error and there effect was outlined by James et al. |^ 
which was derived for the qubits but is equally valid for 
qudits by simple substitution and appropriate change in 
the summation ranges. 



IV. QUTRITS 

As a specific example of how we can implement higher 
dimension tomography consider a qutrit, d — S dimen- 
sional, state. We can write this as 
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where the Xj are now the SU(3) generators and an Iden- 
tity operator Aq. For SU(3) the set of generators are 
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(24) 
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which have been determined using the definitions of 
Eq . ( [l6[]lq ) and the corresponding elementary matrices 

ofE^:fi|j- 

Once we have the expectation values for these opera- 
tors then the density matrix can be reconstructed in the 
same way that it was done for the qubit in Eq.(0): 



P3 = 



|((Ai) 
i((A4) 



+ V3(A3)) 

*(A2)) 
*(A5» 



((A4>-*(A5)) 



^((As) - V3(A3)) |i{X,)-^{X,)) 
l{{Xe)+i{Xj)) 1-V3(A8) 



(25) 



r 



The most direct way to do this is to measure the ex- 
pectation values for the A-operators. However if this is 
not possible let us assume we can measure some set of 
basis states. Consider an arbitrary, but complete, set of 
basis states with the associated projection opera- 

tors {fli — {ipi) {ipil} . These can be linearly related, via a 



cP X (P matrix A, to the A-matrices, 



can thus consider measurement outcomes. 
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E^^^A,. We 



= JVj2AiTT[pX,] 

j=0 



(26) 



where Af is again a constant that will be dependent on 
experimental factors such as detection efficiencies. So we 



find, rj = Af ^ X)f=o(^i ) ^^"^ finally. 
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^ Tli X i 
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In this way the state is reconstructed from the measure- 
ment outcomes in some arbitrary basis and the A-matrix 
which relates the measurement basis to the A-matrices. 
This ^-matrix will be invertible if a complete set of to- 
mographic measurements are made, ie. if we measure in 
a complete basis. The A-matrix becomes the identity in 
the case where we use the generators. 

Take a physical realization of a qutrit in an linear op- 
tics regime. Fi g.(|3| ) shows one way in which a qutrit 
may be realized Jllf . The modes correspond to a photon 
taking the short medium or long paths of the interferom- 
eter. The values of the reflectivities of the beamsplitters 
are such that an even superposition state is generated. 
By varying the phases 0i and 02 a complete basis can be 
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generated, 
|0) 



|0) + l)+2) 
\0) + a\l) + a^\2) 
|0> + |l>+a|2) 
|0) + |l)+a^|2) 



|1) 

\Q) +a^\l) +a\2) (28) 
|0) + a|l) + |2) 
\0)+a^\l) + \2), 



where a — e^^'-l'^ . We can then utihze another 3-arm 




R = 2/3 



FIG. 3: A linear optic implementation of a qutrit where the 
levels of the system are encoded in the pathlength travelled. 
The reflectivities for the beamsplitters, R, are given and the 
extra beamsplitter in the short arm effectively balances the 
superposition of the output. Phase elements in two of the 
arms provide the ability to consider any balanced superposi- 
tion state in the qutrit space. 

interferometer as that shown in Fig.(^) to rotate and 
perform projective measurements on the qutrit. There- 
fore one can perform a series of these projective mea- 
surements and, via the procedure outhned in Eq. (|2^p7|) , 
reconstruct the qutrit. 

The same procedure applies regardless of the architec- 
ture provided we measure a complete set of states. To 
take another optical example, orbital angular momentum 
could be used to realize qutrits (and indeed, qudits), with 
holographic plates generating the qutrit superpositions 
and holographic interferometers acting as analyzers. 

If we now further extend this to two qutrits, which may 
be entangled. 
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we can consider operators of the form \j ( 
related operators. 



(29) 



) Afc , or linearly 



k,l=0 



AfjAfc^Ai, 



(30) 



where the i,j label the rows and fc, I the columns of the 
A-matrix. There will now be cP" — 1 = 80 measurements 
to be made. Therefore, as we did for one qutrit, we 
can again consider the measurement outcomes for states 



of the form {\ipi) ^ \ipj) = \ipij)} with the associated 
projection operators {/iy = /ti (Xi /tj = |}- 

nij ^ Af{ipij\p\Aj) = JVTr[p(pi^ fij)] 
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k,l=0 



(31) 



fej=0 



So we find, rki = Af ^ i^oi^ij) ' ^^'^ finally. 



P23=AA-i J2 (4')"'^.jAfe®A, (32) 
We can then reconstruct the density matrix for the state 
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FIG. 4: A state generation, or measurement, scheme for qu- 
dits using linear optical elements. The beamsplitter reflectiv- 
ities and phases are in complete analogy to the description 
in Fig.(^) for the qutrit generation. For generation and mea- 
surement the number of elements scale as + 3d. 

using the experimental measurement outcomes, riij, and 
this A matrix. Once we have the density matrix for the 
entangled qutrit state we can then consider questions of 
purity and entanglement. We refer the reader to Q which 
gives a thorough exposition with respect to characterizing 
entangled qutrits that is of relevance to both pure and 
mixed states. 

This change of basis is completely general and allows us 
to consider the reconstruction of any discrete system. We 
can now use: the generators; any orthonormal physical 
basis set; or, more importantly, in the case where we have 
limited access to the state space, a non-orthogonal basis. 

As mentioned previously there is significant motiva- 
tion to study entangled d-dimensional states and with 
the reconstruction of the complete density matrix many 
important state characteristics can be determined. In 
practice however, the dimensions will be restricted due 
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to the complexity in implementing the measurements of 
the c?-dimensional state. 

In the case of generating a qudits using a linear opti- 
cal elements the number of elements required to generate 
and hence also measure these higher dimensional states 
increases rapidly. Fig.(U) shows the general scaling for 
a system to generate qudits in a linear optics regime. 
For this implementation the state generation and mea- 
surement requires -|- 3d elements for each qudit. The 
probability of producing these state scales as (1/2)'^^^ 
and similarly for its measurement. Similar complexity 
issues will be relevant regardless of the architecture. 



multipartite-multidimensional systems with limited ac- 
cess to measurements in the Hilbert space. The specific 
example for the qutrit highlights the similarities and dif- 
ferences in going to higher dimensions whilst construct- 
ing an intuitive framework for the Quantum Information 
experimentalist to work. Primarily it is hoped that we 
have made QST relevant and accessible to a wider cross- 
section of the physics community. QST can provide a 
powerful tool for the experimentalist in QIS regardless 
of physical implementation, be it ion trap, quantum dot, 
flux qubit, or photon, to name but a few. 
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